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On the possible origin of gamma ray bursts as a result of 
interaction of relativistic jets with the soft photon field in 
dense stellar regions 
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Abstract. We examine and develop the model of gamma 
ray bursts origin proposed by Shaviv and Dar (1996), ac- 
cording to which the strong gamma ray emission is pro- 
duced by the interaction of the baryonic relativistic jet 
with the soft photon field in a dense stellar region. By the 
simulations of the burst profiles we show that these pro- 
files are sensitive to the jet geometry. Also we derive and 
discuss the model restrictions based on the event rate. 
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1. Introduction 

The idea of cosmological origin of gamma ray bursts (Pri- 
lutski & Usov 1975; Usov & Chibisov 1975) is weU sup- 
ported by their observed isotropy on the sky and non- 
uniform spatial distribution Meegan et al. 1992). The 
enormous luminosity needed for explanation of gamma 
ray bursts if their origin is cosmological imply that the 
best candidates for their sources are the mergers of bi- 
nary neutron stars and/or neutron stars with black holes 
at redshifts (2 ~ 1 — 2), (Blinnikov, et al. 1984; Paczynski 
1991, 1992). Lipunov et al (1995) showed that the cosmo- 
logical origin of gamma ray bursts is consistent with the 
observed B ATSE gamma ray bursts log TV-log S distribu- 
tion, but the beaming of the radiation is required for the 
balancing of the observed and predicted event rates. 

The cosmological models for gamma ray bursts have 
not yet been proven; moreover, they come across severe 
problems, one of them being the problem of the effi- 
ciency of transforming the gravitational binding energy 
into gamma-rays. The explanation of the observed pro- 
files of the bursts remain another unsolved problem. 

Recently Shaviv and Dar (1995, 1996a, b) proposed a 
model of the gamma ray bursts origin which was the first 
one to reproduce the observed burst profiles, and also 
solved the problem of energy transformation into gamma 



rays. In this model, a gamma ray burst is produced by ree- 
mission of the soft field photon {E ^ 1 eV) by the atoms 
of a relativistic jet (or expanding shell) (with F ~ 1000), 
which transforms these photons to those of the energies 
- T^E - 1 MeV. 

It is suggested that such a jet or shell can be generated 
during a binary neutron star merger. In a dense stellar 
region such as a globular cluster core or a galaxy center, 
the particles of the jet would pass close to several stars. 
The light photons filling the vicinity of the stars in the 
comoving frame of the jet would have energy TE ^ 
1 keV. Shaviv and Dar show that these photons could be 
absorbed by photoionization or photoexcitation of heavy 
atoms (like iron) that may be present in the jet, the cross 
section of such absorption being much higher than that 
of Compton scattering. When reemitted, in the rest frame 
of the observer the photons will have energy of the order 
of T^E ~ 1 Mev and will be beamed into the l/P angle. 
Thus, when passing by a star, the jet would produce a 
burst of 7 radiation. 

These 7 photons generated closer to the jet source 
should come to the observer earlier than those gener- 
ated far from it because the jet expansion velocity v = 
C\J\ — l/F^ is smaller than the velocity of light. There- 
fore, each of the stars on the way of the jet would produce 
a peak in the burst profile. Below we will determine the 
shape of this peak. In general, the observed burst profile 
would map the distribution of the soft photon density on 
the way of the jet: u{TC) — ^ F(t), where R is the distance 
from the source to the current point, t is the total time 
taken the signal to pass from the jet source to the ob- 
server: i — R/v+ {D ~ R)/c, where D is the distance from 
the source to the observer. 

The above model explains many of the observed fea- 
tures of gamma ray bursts, including the profiles. 
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2. The geometry of the emitting region 

Let us now discuss the formation of the burst profile with 
taking into account the 7 emission generated out of the 
line of sight. 

If the jet particles moving with the relativistic velocity 
emit isotropically, then in the observer frame their emis- 
sion will be concentrated into a narrow beam of the width 
^0 1/r centered with the particle motion direction. For 
r ^ 1 the angular distribution of the emission in the ob- 
server frame is approximated by 



A{e) = i/{i 



cos^ ( 



In the most general case we can determine the burst 
flux profile by computing the following integral: 



F{t)= J p{R,t-^-^—^)u{R)A{9)dR, 



(2) 



where p{R, t) is the density of rate of the soft to hard pho- 
ton conversion by the jet atoms in space and time, u{R) is 
the density of soft photons and A{6{R,t)) is the angular 
distribution of the emitted light in the observer frame, the 
integral being taken over the whole space. 9{R,t) here is 
the angle between the particle velocity direction and the 
direction to the observer. 

If we assume that the jet is optically thick (for the 
photons of the energy TE in its rest frame) , then the we 
would see the emission from a geometrically thin region 
(the jet "front"). Then we can approximate p{R,t) as 



time. This effect is well known as the rings around super- 
novae. The 7 photon arrival time from the point with the 
coordinates (_R, 9) for the spherical front geometry is 



R (D-Rcos( 
At^ — + 



(7) 



We detect the 7 photons from the surface with At — const 
simultaneously. For the plane jet front, this surface coin- 
cides with the front while for the case of the spherical front 
this surface is an ellipsoid with the jet expansion source 
in a focus: 



(1) R 



vt - DB 
1 - B cos 6* ' 



(8) 



This ellipsoid is strongly elongated in the direction of 
the observer - its eccentricity is as close to unity as 
B « 0.999995. Due to the fact that the jet 7-emission 
is concentrated in a narrow beam with width 0o ~ 1/r, 
we will see the emission only from the part of the ellip- 
soid located inside the l/F cone with the vertex at the jet 
source. It is easy to show that the ellipsoid and the cone 
intersect at Rmin ~ 0-^Rmax, where Rmax is the apocen- 
ter distance of the ellipsoid (Fig. |^). 
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p{R,t) = diR-R{t)), 



(3) 
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where R(t) defines expansion of the jet front. Shaviv and 
Dar assumed that the jet has a plane front moving in 
the direction of the observer. Then in equation (^) 6* = 
anywhere and 



R{t) = vt + R±, V = const. 
Then the observed flux would be 

F{t)^u{{--t)- , 

c 1 — B 



(4) 



(5) 



where B = v/c, i.e. the burst profile reflects the soft pho- 
ton density in the vicinity of the jet. 

Let us now study another case, when the jet front ex- 
pands spherically and then discuss the difference between 
these two geometries of the jet expansion. In our case 



R{t) = vt, V = vn. 



(6) 



The main difference between the spherical and plane 
front is that for the spherical one the emission from dif- 
ferent parts of the front arrive to the observer at different 



Fig. 1. A sketch of the model: We see simultaneously the 7 
photons coming from the ellipsoid expanding after the spherical 
jet front. 



Thus we detect simultaneously the emission coming 
from rather big region of space, which definitely leads to a 
kind of averaging of the photon distribution u(R) details. 
So the details are not so sharp as in the plane front case. 
Their characteristic width should be determined by 0.5i?, 
and not by the distance of the star from the line of sight. 
Calculations show that only the stars located inside or 
close to the l/F cone can produce sharp narrow peaks 
(Fig.|. 

It is important that in the case of spherical jet ge- 
ometry the peak is asymmetric - its tail is much longer 
that the front - which hardly corresponds to the observer 
profiles of gamma ray bursts. 

In order to make the computed profiles consistent with 
the observed ones, the average distance between stars in 
the vicinity of the jet should be 



AR' 



1 -B 



-At 



(9) 
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Fig. 2. The sample burst profile produced by a single star at different angular distances from the explosion-observer line: 
^ = 0, 0.5/r, 1/r, 2/r. The flux is shown in arbitrary units; time is in the units of R{1 — B)/c. 



where At is the average time between the peaks in the 
burst. For T = 1000 AR - 0.02Ai pc/sec. It means that 
for reahstic gamma ray burst profiles the star density 
should exceed 

~ AR-^ ~ 10^ pc"^ (10) 

Such high density can occur only in globular clusters and 
in the galactic nuclei. The size of the region where the 
burst is formed should be of the order 1 pc which is similar 
to that of the globular cluster core. 

After the above the gamma ray burst profiles can be 
simulated from the integrating of the contributions of all 
the stars in the vicinity of the jet. We assumed that in the 
globular cluster core 10^ stars are located within 1 pc^. 

3. The event rate 

In order to bind the gamma ray bursts to any events at cos- 
mological distances the rate density of such events should 
be 7?. ~ 10~^ Mpc~^ yr~^. If we allow for a beamed radi- 
ation, then TZ ^ 2 - 10~3/3^ Mpc~^ yr-\ where /3o is the 
beam width in degrees. 

The rate of neutron star mergers in the dense stellar 
regions - globular clusters and galactic nuclei - is still not 
clear. 

Numerical simulations of the evolution of close binaries 
(Lipunov et al, 1996) show that the neutron stars merger 



rate is strongly dependent on the stellar population age. 

For the old stars of the globular clusters it should not 
exceed lO^'"" yr^^ per IO^^Mq; for the relatively young 
stars in the galactic nuclei it should be of the order of 
10"'' yr~^ per 1O^^M0. After taking into account the fact 
that the mass fraction of these dense star configurations in 
the galaxies is not more than 10^^ wc obtain the optimistic 
upper limit of about r = 10^*' yr^^ per W^^Mq galaxy. 

For the flat Universe with ilgtars = 0.005 and Hq = 
75 km/c/Mpc it corresponds to the rate density TZ = 
10~^ yr — 1. So, there is no place for the beaming, and, 
thus, for the jets. At least, the jets should be very wide 
(which supports the spherical geometry case described in 
the previous section, in fact, if the jet front is plane, all 
emission should be in very narrow l/F beam). 

If only the neutron star merger rate in globular clusters 
or galactic nuclei is essentially higher than 10~^ yr"'^ per 
IO^^Mq the above model can be used for explanation of 
the gamma ray bursts origin. 

4. Conclusions 

We have studied some of the conditions necessary for the 
explanation of gamma ray bursts as the interaction of the 
baryonic relativistic jet with the soft photon field in the 
dense star regions (Shaviv, Dar 1996). The main attrac- 
tiveness of this model is that it is the only one that pre- 
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Fig. 3. The example simulated gamma ray bursts profiles computed with a real distribution of stars in a globular cluster. 



diets the observable gamma ray burst profiles. We have 
found out two difficulties of these model, none of them 
being fatal but each of them putting strong restrictions 
on it. First of all, the main achievement of the model - 
the reproduction of the observational burst profiles - has 
seemed to be strongly sensitive to the jet front geome- 
try. The shape of the peaks in the spherical case becomes 
more wide and asymmetric, and the number of the peaks 
becomes smaller. 

The second is that it is hard to agree the theoretical 
event rate with the observed one. The neutron star merger 
rate in globular clusters and/or galactic nuclei should be 
high and the beaming of the jet should be small. 

The spherical geometry is more consistent with the 
event rate than the plane one because it allows for beam- 
ing wider than l/F, but is in worse agreement with the 
observed bursts profiles. The detailed investigation of the 
jet expansion geometry would allow now for the further 
treatment of this model of gamma ray bursts origin. 
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